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■ ....summary / ■ :, ' : . 

•As the initial. step in the. analysis of stress ,di.Sr. 
tribution in three -dimens i onally curved rings (as employed 
as stiff eners in. stressed skin aircraft ..des igns ) the. ring 
(fig., l) formed "by the intersection of tto circular cyl- 
inders is explored for three categories of. load: tension 
in both cylinders (produced "by hydrostatic pressure, on the 
cylinder walls), axial force in the large cylinder, and 
lastly, shear in. the large cylinder. The discussion of 
these three load. cases enables general conclusions con- 
cerning the behavior of the ring stressed by • the .• she 11 
forces and affords numerical data for the mos t . important 
load categories (obtainable from the computed by super- 
position). The quantitative r e s ul-t s. ar e". i 1 lus t rat ed in 
figures 12, 14, anri 15, and condensed .in simple approxi- 
mate f ormulas _ through (4.9), (6.. 2), and • -('7 . $ ). . Qualita- 
tively, it can be- stated that , on wings which do' not de- 
part excessively from the plane,, the moment M 2 about 
the normal axis (hence that of the three moments which is 
other than zero even on a perfectly straight wing) remains 
the paramount stress; and not until there is. a very ap- 
preciable three-dimensional curvature (when the ratio 
b/a of the- cylinder radii approaches 1) do the two other 
components of the three-dimensional moment vector, the 
bending moment M 2 and the torque Mj .become perceptible. 
Since Mg , as the. graphs indicate, varies but little 
with e = b/a if suitable reference quantities are chosen, 
rings for -"small" openings can be computed as straight- 
rings with very good approximation. 

The (closed) ring is 'statically indeterminate. It 
effectively evades an excessive stress induced by insuffi- 
cient torsional stiffness- hy responding to the load largely 
with bending moments M s — still, on rings fully ineffec- 
tive in torsion, it is recommended that the existence of 
the shear stresses 'within permissible limits be confirmed 
•by approximating with the help of the cited empirical formu- 
las'.' ■ 

'-*" Spannungen in Aus s chn i 1 1 ve r s t e if unge n , " Luf tf ahr tf or s chung , 
vol. 18, no. 7, July 19, 1941, pp. 253-61. 
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The three explored load dis tr ibut i ons - are three 
column loads for the unstiffened cylindrical shell; hence 
they create in the undisturbed shell a pure membrane 
stress attitude. The calculation -is predicated on the 
assumption that this membrane s t r e s s .at t itude is net ma- 
terially disturbed by the elastic interference effect be- 
tween the stiffener opening and the skin. This assumption 
is met in the "extreme" case of a very stiff ring and a 
thin-wall shell without frames (or with frames located at 
some distance from the opening). 3Tor the frameless shell 
would, have to attempt to terminate the "interference loads" 
returned by the ring through cross stresses and bending 
moments; since these do not become large in the thin shell 
and are damped quickly besides, any 'aid" of the shell for 
the ring can manifest itself merely in the formation of a 
small effective border zone which takes nothing essential 
away from the ring. 

In the opposite extreme case (not discussed here) of 
a ring rigid in' strain but flexible in bending and of a 
shell closed all around by closely placed stiff frames or 
curved floors — "egg surface" - the state cf stress and 
strain is. utterly different. Shell and ring are for the 
most part subject to diaphragm and axial stresses, and 
stressed in bending solely by the constrained stresses due 
to the. incompatibility of the form changes. : The case is 
of little practical concern, since structural reasons usu- 
ally .call for rings which are far from inef f ect ive -in bend- 
ing. 

The. true' shell lies between the' two extremes. If the 
ring is distinctly rigid in bending and the shell either 
is thick— walled or forms an egg surface, a complicated elas- 
tic interference' effect results which defies calculation' 
and must be ascertained experimentally. The present 
solution supplies the basis for such experiments by enabling 
the estimation of the "maximum bending stresses to be expected 
through -'the determination of their upper limit. 



INTRODUCTION 
The Flat Ring 



Concerning the exact stress distribution of three 
dimens ionally curved r ings , such as are used as stiffeners 
on openings in shells of all kinds, little data are avail—, 
able. The present study treats as a typical example a 
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ring the center.; line -.of . which is' produced "by the 'intern- 
section of two circular cylinders'" of. different- diameter.:* 
Three load cases are- analyzed.: 

1. Axial and circumferential stresses in both cyl- 

inders, ,t;he cylinder /stresses themselves tr: 
he in the ratio conformal to the cylinders 
loaded under internal pressure 

2. Pure longitudinal tension in the- large cylinder 

3. Pure shear (torsion) in the large cylinder 

To simplify the calculation, it is assumed that the 
ring, compared to the shell, is very strong, so that its 
deformations have- no perceptible effect on the stress 
c ondit i on- in the shell. -.This provides an. upper limit for 
the' ring stresses actually produced in a shell design, 
for, according to the theory of stressed skin statics the 
shells, by elastic flexibility of the ring, regroup the 
forces deposited on it in such a manner that the ring is 
relieved. 

I±°!t^_£sLi=!L_l • ~ The solution can be given immediately 
in the extreme case a»b (figs. 1 and 7), that is, 
if the ring is "practically" flat. Then the forces ex- 
erted by the small cylinder are secondary alongside those 
of the large cylinder: the force pa along the circum- 
ferential circles and the force pa/2 along the generat- 
ing axis. Since an equal tension pa/2 from all sides 
stretches the ring without twisting it, the bending stress 
can be computed as if the ring were loaded in the manner 
shown in figure 2. (Load cases 1 and 2 become identical 
except for the exchange of axes.) 

The equilibrium conditions on the ring element ds = 
bdcp are expressed in vectorial form from the' very 'start 
in view of their subsequent application to the three- 
dimensional problem,. Thus t "denotes, the unit vector of 
the tangent pointing toward increasing arc length s, n 
the unit vector of the normal toward the center of the 
circle, k the unit vector at right angles to the- plane 
of the ring (fig. 3); 1 is to indicate the resultant, 



*In fig. 1 the smaller cylinder -is shown outside the large 
one. But the ring formulas apply exactly, if extending 
wholly or partly in th'e large cylinder; the small cylinder 
can be arbitrarily s.hprt; it can be formed by the ring 
itself, for instance. 
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M the moment of the section stresses applied at a sect'isn 
with the outs ide normal t; at the section -"boundary for 
which t is the inside normal the resultants — N , -M are 
effective. Since £ and M vary with s, the amounts 
on "front" and "rear" of a piece of length ds differ "by 
d¥, dM/ In consequence, the force equi libr iunr spec if ies 
according to figure 4: 

dl + £ ds = 0 

the moment equilibrium (in absence of external moment 
loading) 

dM + (t ds) x I = 0 

(the choice of moment reference point within length ds 
"being immaterial, since t he . d if f e r enc e s of higher order 
accruing therefrom "become small within the limit d— >-0). 
These vectorial equilibrium expressions for the bar element 



— ■ + p = 0 
ds ~ 

dM 

— = + t X N = 0 
ds 



(1.1) 



can, if load, sectional force, and moment are divided into 
components along Jb , n , k , be written in the form 



(N t. + Q a) + p t t + p n a = 0 



d 

ds 

= d - (M k) + t' x (N t + Q a) = 0 
As _ _ 



and because of 

dt 1 dn 



n , £S = _±t,tXn=k (1.2) 
ds b ds . b . . 



are equivalent to the scalor 

dW Q dQ N : " dM 



= - Pt 



— + 



ds b ds b n " ds 



(1.3) 



which, in this two-dimensional case could naturally have 
been read off as" well fr.om figure 5. - .. . 

In this particular load study the components p^ , p n 
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(dimensions: force per unit length) should be replaced by 

p t =s ~ sin cp cos cp, ' -p n = sin 2 9 (1.3 ') 

according tofigure 6, The integration of (1.3) presents 
no difficulty. With ds = bd9, the first two equations 
give: 



N = A cos cp + B sin 9 + 2~ c 0 s 2 cp 
Q, = B cos cp — A sin cp — — ~ cos cp sin cpj 



(1.4) 



The two integration constants A and B follow from the 
symmetry requirements according to which the cross stress 
cp = 0 and tt/3 must disappear. Then 

A =. B = 0 (1.4') 

Entering (1.4) and ("1 . 4 1 ) in (1.3) and integrating affords 

M = X — cos 2 cp (1.5) 

Integration constant X remains statically indeterminant ; 
it can "be computed "by means of Cas t igliano 1 s principle of 
least— strain energy 



dX J 



M 2 
2 E J 



ds = 0 (1.6) 



For the specific case EJ = const, equations (1-5) and 
(l . 6 ) give ' ' 

M = — cob 2cp (1.7) 

8 

In consequence , the two extreme values' of the moment at 
•points cp = 0 , cp = tt/2 are inversely equivalent and 
amount to 

2 



pab ' 
T 



M max = ~ — (I-?') 



Equation (l.v5) for the moment can equally be derived 
by another process which is much more simple-, to wit: 
'According to figure 6, the first vector equation (l.l), 
when resolved along the vectors i_, ,j, k (place inde- 
pendent) charac t er is t iz ing the space directions x , y, z, 
instead of along the "natural" variable directions t_ , n, 

k read 
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dN pa 

— . 1- ___ j s in cp = 0 

tdcp 2 1 

the integration of this equation is even simpler than that 
of (1.3); -we. get 

pah pah 
I = 5 0 + -g- i cos <? = Cji + C 2 j. + — j cos cp (1-8) 

that is , 

pah "[ 
Q = N x n = — . Cj cos cp — C s sin cp - — -- sin cp cos cp 

pah / (LB' 

N = N x t = - C j sin cp + C 2 cos cp + - — c o s 2 cp 

2 j 

Integration constants G x and C 2 disaunear because 
0.(0) =Q, (tt/2) = 0; hence 

pah 

I = — j cos tp (1.8") 

and with it follows, because t x j = — k sin cp and Ms M k 
from the second equation of (l.l), and equation (1.5) 

pah f pah s 

M = — ~ / cos cp sin cp d cp = X — cos' cp (1.9) 

2 J 4 

The calculation of the three dimens i cnally curved 
ring also proceeds in two stages: the determination of 
the intersection resultants (forces and moments) as far as 
is possible on the basis of the static statements (l.l), 
and the solution of the statically indeterminate quantities 
on the basis of the strain conditions. The problem is 
most easily solved if the vectors for the integration of 
the differential equations (l.l) are resolved along the 
place— independent system of unit vectors i , j. , k, and 
then the transfer to a system of axes attached to the 
space curve (tangent t_ and two normals n 1 , n 2 carried 
out the prediction of the integration constants from the 
strain- and symmetry conditions'— and these only — necessi- 
tates resolution along the natural axes ;t , n lt n 2 ,' be- 
cause the strain law and the symmetry expressions are 
amenable to simple formulation only for such components of 
the force and moment vectors. 



UACA- Technical' Me'mora'ndu-m" So. -1-005 



7 



GEOMETRY OE THE . SPACE CiJEVE = 

The first step in solving the three-dimensional ring 
problem is the lay ing_ down „ pf the f ormulass. character i z ing 
the geometry of t he" spa-ce curve,, — With the notation of 
figure 7 the space curve is given by the two' formulas 



*2 + y2 = bs, y2 + z E = a?. 

With the ch,o:ic:e:. of .-the angle -<p pr o je ct e'd Int o •. the xy 
plane as place parameter and b/a = the triple: equation 

reads 



x y . z 

— = cos cp , — = s m cp , — 

b b b 



sin 



9 



(3,1) 



'. Erom the geometry' of the' space curve" represented ' by . 
(2.1) two groups; of formulas ' are applied: .1) the expres- 
sions for the arc length and for the three unit vectors 
of an "accompanying triangle," 2) the relations expressing 
the position of the curve element with respect to the 
force directions i — The arc length follows from 



ds = 



f d s N \ 

V dcp y 



d 9 



/ dx"\ s '. /dy^ /dz\ 2 



\dcpy 



+ • + 

vdcp/' : 



dcp 



with 
dx 
dcp 



dz 
dcp 



at 



— X = — b sin cp , — - y = b cos cp 

.dcp • •. 



- = i = ^ -b 



sin cp cos cp_ 



JTT t . 



3 sin 3 cp 



ds = bdcp 




• " S '. 4 ■ 

— e sin cp 

v •' 2 r . ' 2 ' T 
1 — € s m cp 



(2.2) 



(2.3) 



The introduction of " the "angle \J> ' of - the -space 
ge'nt" with respect' to' -.it's projection ■ 



curve tan— 
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dx dy dz 

— i + — j + •— k 
ds ds ds 



- <kl +. 73 + *Jk)s— = Ui.+ yj. + z&) 



ds ~ ■ 15 

with cos and sin of \j/ "being given "by 



, 1 - e s sin 2 cp € sin cp cos cp 

cos- \|/ = / ■ -, Bin ♦ = — -~— ' (2.4') 

1 — e 8 sin- cp vl — £ sin cp 

according to (2. 4), 

Defining the two mormals n x and n 2 "by the stipu- 
lation (of itself arbitrary, "but in view, of- the simplicity 

of .the formulas appropriate), to place n a "horizontal;'! 

(n. 1 xk=D, n.iXi=0, n 1 =l),. the formulas for the three axes 
read* 



t_ = (— i sin cp + _j ces cp) cos- \|/ .+ k sin \|/ 
.'n 1 = -(j.cescp + _jsincp) 
t X n x = n 3 = (i sin cp — j cos cp)sin \|/ + k cos \|/ 



^ (2,5) 



The forces that stress the ring are applied on it "by 
the cylinder skins. 

Consider figure 9, which represents a ring element ds 
and (slightly shifted) a skin element at one side of which 
(with the normal v) a shear force. K is applied. The 
load e acting on the element of the ring is , "because- of 
the equilibrium in the skin element, given by 

p = K cos 'a - K(a"v) (2,6) 

There the ring normal n is characterized by the fact that 



*The use of the so— called natural axes t_ , n, b (tangent, 
principal normal, binormal) for describing the curve is 
unnecessary and, in general inappropriate.,' For' a deter- 
mination of the natural axes which has not implicit connec- 
tion with the principal inertia axes of the section requires 
the' knowledge of the third ■ derivat ions of the- system (2.1); 
whereas (2,5) follows from the first d'erivations only, The 
fact that one of the curvature components disappears on the 
natural axes is an advantage which is of no consequence com- 
pared to this drawback. 
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it with V and t is located 'in one plane; hence it 
falls in the tangential plane of the particular cylinder. 
It is most simply obtained over the surface normal Hj 
or nj.j ., • respectively, to which it must he at right an—. 

gles. -The outside normal of a circular .cylinder falls 
along the radius vectsr. Therefore, "because of (2.1) 



= j_ e sin cp + k 



J 1 - c 



1 



s in° cp 



.= i cos cp +. j sin cp.( = — n 1 ) 



(2.7) 



With the identifying signs of figure 10, the desired curve 
normals "become 



n a = i X- n T = i 



c os cp 



sin cp 



+ j gf-^lsljglglS^g - k £ sin 2 cp cos *j 



(2.8) 



as- affecting the large cylinder 

n- b (= n 2 ) = njj X t^ = i_ sin cp sin \|/ 

— j 'cos ip sin i' + k cos ^ ' 

relative' to' the small' cylinder. 

.Later on, the angle "between the tangent' 

_t' 1 = (j yi — €2 sin 2 cp — 3c 6 sin cp.) 

at a circumferential circle of the large cylinder and the 
normal n & - is particularly needed; • 



cos(jb',n a ) = t n a = sin cp cos\l/ 



(2.8 ») 



THE EQUILIBRIUM' EQUATIONS 



The' eq.uilibr ium equations for the ring element (dis- 
tributed outside moments discounted for the time being) 
read in vector form as in -the two-dimensional case: 
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dH dH . , 

'M + £= 0,-=+t_x.N=0 (3.1) 
ds "" ds 

with _pd s as vector of the external force applied at the 
ring element, as results frotc (2.6); the shear force vec- 
tor - 

M = Nt + Q, 1 n l + ^ 2 n 2 (3 . 1 1 ) 

has this time the longitudinal force N and the two cross 
forces Q, 1 and Q, 2 as components ; ■ the shear moment vec — 
t or 

M = Mjt + Mjnj + ;! 3 n 2 (3.1") 

has the torsion moment and tending moments K 1 and 

M a , ~ ■ 

Equation (3.1) is integrated in two stages: 

r s -co 

1. N » H(0) - / p ds = 1(C) - 0/ p_ -— -r (3.2) 

/ ^ « 0 cos V 

. m = M(0) - f S t x Eds = M(0) - b / US ~~- (3.3 ) 
4) 



J. " ~ v' 0 CCS V 



The integration constant IT ( 0 ) is again found by sym- 
metry considerations, one component of the second constant 
M(0) remains indeterminate. The symmetry of system and 
load requires the disappearance of the three ant isymmetr i — 
cal quantities Q, x , Q, 2 , and (the shear resultants) at 

points cp=0 .and cp=tr/2;from l ^ 1 (0) = 0 and 0,^(71/2) = 
0 follow the i and j_ components, from Q, 3 (0) = 0,3(^/2) 
0, the k component of U(0); from Mj (0) = 0 and 

%(£) = 0 the i and j components of M(0) - the k 

component of M(0) — the k component of M(0) remains 
to be determined by a strain (equation.* 

Denoting the intensities of the shear load of the two 
cylinders along the generating axis and the circumferen- 
tial circle with 



*Since, for reasons of symmetry, N and M contain only 
the even— number harmonics in cp , Q, 2 (rr/2) = 0 follows from 
Q, 2 (0) = 0. Hence the cited 6 symmetry condit-ions yield 
only 5 independent equations for the determination of the 
6 integration constants. 
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Pi • -Pa for the large cylinder 
P3 » P4 f ov the small cylinder 
equations (2.6), (2.8), and (2.8') afford 

. COS (B . j -. 

P- = Pii ------------ + PsQ VI - e 2 sin 2 cp 

VI - € 2 sin 4 cp 



6 k sin cp) 



Xsin cp cos \[/ + p 3 k cos \J/ + p 4 (j cos cp - j. sin cp) sin ^ 

(3 ..5) 

whence 

- [ Pl i- 7== =~L: = 0 3 = - p^i / __ <!££-!£_ 



*1 



Vi - £ 2 sin 4 cp 



dcp 



J 



£ 2 sin 2 cp 
= — PihJ. — arc sin( sin cp), 



P2J 



Vl~ 



£ 2 sin 2 cp sin cp cos\l' ds = 



+ — - — j coscoVl — e sin° cp + — (l _ e "*) sinh 



™ /l <r 2 • 2 , 1 /•, 2\ , , - , C COS If 



AT 



?38 a y 








s«- - 





£ sin 2 cp cos v|/ ds = p 2 k e b 



2 _ s in cp cos cp 
2 2 



£ sin 2 cp cos cp 



£ 2 sin 2 cp 



dcp 



= ~P 4 i ~ [ si * cp^ - sin 2 cp - I (sin' 1 ) ( € sin cp) 1 • 



p 4 j, cos cp sm 



,1; _ ^ v • P sin cp cos 8 cp 
W as = p 4 l)j e / -y======t : == =: = === x == a cp 

J VI — £ 2 s in 2 cp 



r - ^ _ 

P 4 1 ~Uos cp VI- e 2 sin 2 cp-i(i- e 2 )sinh~ 1 4==^== 
2eL VI - ;J 



e 2 - 



(3.6) 
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It is noted that the load portions p 2 and p 3 , . 
which produce N 22 and N3 , of themselves for? no equi- 
librium groups , because the N contain a non— periodic 
port i on .which cancels out when p 3 is equated to 
- p 2 .. This occurs, for instance, if the four forces 

originate through the same internal pressure p in the 
two cylinders., 

pa pb 
Pi = --. P 2 = Pa, P 3 = ~. P4 = P^ 

Those values, written in (3.6) and condensed, give 



(3,7) 



e 2 s in 2 cp 



N = i gin cp Jl 

+ 2j cos cp 1 - e 2 s in 2 cp - € k sincpcoscp] (3,6) 

Force N is represented by (3.6') with the correct in- 
tegration constants, because the i and k portion dis- 
appears at cp = 0 , the j portion^.at <p = tt/2 . 

The solution of the moments according to (3,3) is 
predicated on the three vectors t x i , t. x j_ > x k. 
According to (2.5) 



t x i = _j sin"v(j - k cos cp cos , 

t x i = -i sin ^ -k sin 9 cos' it/, } 

i 

t_ x k = i_ cos cp cos \|/ + i. sin cp cos 



(3.8) 



which, entered along with (3.6') in (3.3) and integrated, 
gives with M(0) = Xk 

p a b 3 / e cos cp (1 e sin cp) 

M a X k + — i , - k — 

2 V~ 3 3 £2 



(3.9) 



By means of the transformation equations (2..5) the natural 
components Mq> , M x M 2 of the moment vector then follow at 



M T . = _ 



Mj = 



pab 2 1 - e 2 



6 e 
pab 3 



sin cp cos cp cos\J/ + X sin \|/ -, 



£ cos cp , 



> (3.9') 



pab 3 (1 - e 2 sin 3 cp) s + e 4 sin 3 cp cos cp 
M a = - r .~ — . + X cosv|/ 



A' 



Sin cp 
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these expressions satisfy, as is seen, t he., s ymme t ry con- 
ditions M T (0) = M t (tt/2) = 0. 

The course of the three moments, particularly in the 
important practical case of c <s; 1 (small openings ) , is 
of interest. Expansion in powers, of "e affords 



pah 2 1 3 
M T = — - — — sin cp cos cp | 1 — e 



2 

— e X sincpcos cp( 1 + cos. cp + 



/ , s in 2 cp c os 2 cp \ 1 

C i + — v-— ;-••] 

•)■ 



2 

M x = - £ C os 4 cp. (3.9") 

pab 2 r ~ / 1 4 \ P. 

M 3 = — ^- ^ -B- - (2 s in 2 cp - - s in cp^ . . V 

. v ^-i £ 2 • 2 2 ^ 

+ X 1 — — sin cpcos cp— .... 
V 2 J 

SOLUTION OF THE STATICALLY INDETERMINATE X 
Simple Formulas for Maximum Moments 

The prediction of the integration constant X is pred- 
icated upon a strain equation. . If expressed in the form of 
Cas t igl iano * s principle of least strain energy, the geometry 
of the strain condition is secondary (reference 1). In 
vectorial form Cas t i gl ian o 1 s requirement ■ reads 

J~ M k ds = Min (4.1) 

with k the vector of the curvature change 

k^4t_ + K 1 n 1 +K 2 n 3 

with three components: twist and curvature changes- 

k i , K 3 . The thin, slightly curved "bar serves as hasic 
strain law, the general case of diagonal bending being 
analyzed at once. Taking into consideration 

a x = E(w M z + v"y) = E ( i^z K s f) 

(y , z distances from the centroidal fiber of the bar) - 



14 



NACA Technical Memorandum Noi' 1005 



,-.M 

the law reads 



z d F , M ; 



:- - / 



a x y • d T 



M 3 = G J 8 



E J 1S k 2 



12 a 3 



(4.2) 



with J lt J 2 , J 13 the inertia and the centrifugal moments 
referred to axes n_j and n 2 (y and z") , the tor- 

sional resistance. The solution of (4.2) inserted in 
(4.1) gives, with the abbr ev iat i on 

A = 1 r 

J l J 2 



the equation for X: 



d X J J . EA \ ■ J s J x J s / 

r m 
-f -~ 



~— sin 1 -!/ ds + X j 



b / M 3 



~e~ j~~a~ + y 



c os ^ d s 
sin 2 ii ds 



0 J 



T 



] 



eaV 'IT + ^Tjfj d ' + / 13 tan * d * 



+ x 



b cos ^ 4 9' s in 3 ^ 

E~~J 2 _ A + ; 



where .M 2 



M 



(o) 



J COS ^ &Jj 

(o) 



G J j 
d cp 



(4,3) 



are the statically deter- 



mined portions in (3.9'). 



The evaluation of (4.3) is predicated on the course 
of the quantities J with s and . cp , respectively, "To 



"begin with., — -I— - — AS will usually "be small, compared with 

M 2 0) , because M^ o) 
M ' ' ' 



(o) 



is .smaller by ijhree e powers than 



2 and factor J-i 2 /Ji itself ' will . be perceptibly smaller 



than 1, So, the second term of the first parenthesis can 
be discounted for the rough calculation of X. The same 
applies to factor A in the denominator, which, even by 
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marked, departure from 1 , does not affect the result for 
X materially , since it acts in' the same sense in the 
numerator and the denominator, If desired', it can he 
alloyed for . in. the determination of the still remaining 
essential parameter 



a. = — — (< 1 ) . cr • = 

E J g E J" g A 



(4.4) 



which indicates the ratio of torsional to flexural stiff- 
ness. Assume average values for J 2 and Jrp to he in- 
dependent of s, the X equation simplifies to 



with 



X = 



A + B a 
C + D a 



rr /a 



' tt/2 



= J M^j° ^ tan\|/ d cp, 



B 



o 

tt/s 



M^ 0 ^ d cp 



s in 2 \j; 
c os \|/ 



d cp 



o 

rr /a 



cos ^ d cp 



(4.5) 



(4.5 1 ) 



For large € values the integrals A...D must "be numeri- 
cally evaluated, for small e integration "by series ex- 
pansion is suggested. 



A = 



pah 



r 

-J 



.2 a / 2/ 

sin cp cos cp'^1 - c _ 



s in 4 cp *N 



J~" 7 



dcp 



- ^- ( 1 - e' 



B = ^ 



6X8 

~ I ^ l-€ 2 ^2 S in 3 cp-|- sin 4 cpy + , .ydcp 



o 5 



pah 2 / 
6 



/ 



(1 



2 S / € '4. 

sin cp cos cp/1 + — (sin 3 cp + sin cp) + . . .J dcp 



= ii(i 

8. V 



i-(l + 5/4) ) 
4 / 



e 2 sih 2 cp cos 2 cp- ; . . )dcp = (l- ™(1- 3/4) 



V 
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Figure 11 illus trates ".the result- of the ' caleulat ion ■ 
of X for a specified axis ratio (e = l/2), .along .with 
the three stat ically - determinate quantities Mr^ 0 ^, u[° \ 

Mg 0 ^ and - for different a values - the curves X'sinl{/ 
and X cos i> , respectively, from which the final moments 
M T and M 2 . must he marked off. 

At e = 1, the ring degenerates to two flat pieces 
(semi-ellipsoids), which meet under a right angle. The 
particular loading (3.7) which is symmetrical to the two 
ellipsoidal planes, stresses then each half of the ring 
in its plane only. In other words, the statically deter- 
minate portion of the torsion moment must disappear (first 
equation of (3.9')) and the statically determinate portions 
of both tending moments must combine to a bending moment 
ab out the normal to the plane of the ellipse. In point 
of fact, the moment vector at £ = 1 : (a = b) -for the ' 
first and fourth quadrants reads, according to (3.9) 

„ ( 0 ) pab 3 3 . . 
M = cos cp(i - k) 



for the second and third quadrants 



(0) . pab 2 3 . 

M . = cos q> ( 1 + k ) 

6 ' 



hence is at right angles to the plane of the ellipse. At 
a = 0, that is, vanishing .torsional stiffness, the stat- • 
ically indeterminate (4.5) disappears, because the other 
half' cannot absorb a bending moment as torsion moment at 
cp = tt/2. Hence m(°)(tt/2) = 0; - at a £ 0 a reciprocal 
restraint occurs which produces torsion and bending moments 
diverging "from m(o), that is, twists the ring half out 
of its plane. (The result of. the calculation for any £ 
and a illustrated in figure 12, indicates that M £ 0 ' 
and M 3 (°) actually disagree at a £ 0.) 

The maximum amounts of the moments in relation to a 
and e are of particular concern. For e < l/2, they 
are readily obtained by means of the series expansions 
-(4.5"). For X there is obtained 

e 2 

-a a (1-0. 8125-e 2 +0.2275 € ) + — (l-C . 8437 e 2 )+ . . . 

i <«••> ■ 

a (1-0, 0625€ 2 -0. 0175 €*)+-- ( 1+0 . 4063 € 2 ) + ■ . . .■ 
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hence" at 'a f 0 



pab B " _1_ 
6 L € 2 



3 : / ■ 1 : A • 

- + e 2 ( 0.1981 ; + 

4, A 16a/ 



(4.6 ') 



at a = 0 (very low t ors i onal • s t if f nes s of ring) " 



pah 2 "15 

— [■75 ~ 4 + °;- 2738 £2 + • 



(4. 6" ) 



(The statically indeterminate destroys, as it should, in 
both cases the strongest term in the expressions (3:9') 
for Mrp and ; M 3 ). .For the maximum values' of bending mo- 
ment M 2l which are located at the symmetry points cp = . 
0, tt/2 ... (cos \|/ = 1), there is obtained 



at a / 0 
M 3 (cp = 0) = 



pab"* ; 1 1 3 s / ^ 1 
_____ _ + _ + e .. i 0.1981- - - — 

"6 l_ £ € 4 \ 16a 



pab . s i \ 

= : _ | r_ £ : o .264 i + 

8 . L •. . \ . • • 12a y 



M 2 (cp . = tt/2 ) = 



^ i ' 1 , 3 -3e 2 1 3 _. . 

— " — I ~ --s + -— + -- — £ -+-e 3 (. . )+.■ 

6 i_ e 2 8 £ 2 4 



= P^ 3 j~! _ 



0.236 + 



12a 



at ■ • a = 0 

M 2 (cp = 0) = - -Zl-\ 1 _ 0.365 e 



(4.7) 



pab" T 5' 
8 L 3 



M s (cp = tt/2) = _-t_i 



8 L 3 



- - 0.168 C S + 




(4.7 ') 



The equations • (4.7) confirm first the "two-dimensional" 



result 



°) t M maxI 



pab 1 



they further indicate that 



M 2 is reduced at cp _ H and likewise at cp = 0 , as long 



as a > 



12 X 0.264 



= 0.315. For smaller a values the sec- 
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ond "bracket of the first expression reverses, signs: the 
maximum amount of the "up" "bending moment (<0) at point 
0 is greater in the three— dimens ional (e £ 0) than in 
the two—dimensional case (e =? 0). That M 2max must be- 
come greater at small a values should not be surprising: 
the ring tries above all to evade a torsional stress; 
hence it can not devote much attention to the reduction 
of the maximum bending moment (bending and torsion are 
interconnected on t he • t hr e e— dime ns i onal curved ring). At 
the limit a =? 0 this tendency of the ring even results 
in a radical departure ..of the M 3 (a,e) curve from the 
others: curve M 3 (e) m _ n ri = n alone, passes at e = 0. 



through — 



24 



pab 



cp=o ,a=o 



instead of the point — 
drops monot onically to — i pab 2 at e 1 



1 

8 



pab 



as a 



J , . and 
f unc— 



Figure 12, where M 3 (e ) 
shows that the 



tion of € 

different a values 

for some distance to the boundary curve 
small torsional stiffness values, while 



5/24 pab 2 is not reached for finite values of 



has been plotted for. 
M 2 (c) curve clings 
M 3 (a = 0) for 
the extreme value 



a. 



The maximum torsion moment (other than for a = 0, 
where M m is on the whole very small) is approximately 

located midway between the symmetry point 0 and tt/2. 

Hence the amount M t (tt/4) is determined as approximation 

for. Mm .For small values e ' the expansion in series 
•'-max 

is again recommended.. 

pab 2 • f i 

M T = e sincpcos'cp-<;-,T 

6 l € " 

4 



. . s in 2 cp c os 2 cp A b, . 
1 + ±J- e (...) 



-, 1 - 3 + 
e 2 4 



(...) + 



V 



1 + -- sin 4 cp + sin 8 cp 



8 



M T (rr/4) = 2l! i(- i + c 2 (± 
6 2 L 2 \32 



5 

.32 

3 - 



0.1984 + 



M4.8) 



Pb 

= - — II - c 
.24 l 



16a 7 

( — 0.842 ) 

V 8a / 



■i 



At a - 0 the next term in the series expansion along e 
likewise disappears for cp = tt/4, leaving: 

' - 3 

M t (tt/4) = I c ( + 0.2738 

■ (a=o) 12 L V 12 8 32 12 8 J- 



\1 



3 2 



= - p b ex 0.0046 » 0 
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The result of this discussion is ' that , '■ a-s " regards 
,M gl the stress analysis with a small safety margin can 
"be carried out according t o ' t he • s imple formula 

>max = -g-P a *° S : . .". . : < 4 - 9& ] 

For the ' maximum value of the other "bending moment M lt 
equation (3.9') affords ., . . . . 

M ima X = ^ = | pab* ^ ^- (4,9b) 

For Mrp a simple approximation, which is practical up 

max 

to e = 0.7 and remains on the safe side near £ = 1, 
is given by the zero point tangent in-figure 12 

M w - h ^ s " it. (4 - 9c) 

'M 1 and are in our particular load case independent 

• of the radius . of the- large cylinder-,' 

EFFECT OF ECCENTRIC STRESS APPLICATION 



The effect of the moments which stress the ring di- 
rect following eccentric load application is secondary 
compared to the moments (3,9) set up by the cross stress 
because of the small lever arms. However, an appraisal 
seems desirable; it can be achieved by means of equation 

(3,5). To determine the order of magnitude of the addi- 
tional moments an approximate assumption is that the 
forces p j and p B apply at a lever arm h x along k, 
and forces p and p* at a lever arm h along gj 

(fig, 13). Then the localized load (3.5) produces the 
following distributed moments m: 



m = - p ih-^i— = =5= ==~=i= s == + p 2 h! i_ sin cp cos"4>Vl -r e 3 sin 2 ep 
vi _ e 2 sin 4 cp 

.+ P3 h s {-* i sin <P '.+ J c°s'' cp) + p 4 k'sin\|/ (5,1) 

which give rise to the shear moments 

M m = / m ds = b / m — . — — (5.2) 
- J ~* J cos \!/ 



20 



HACA Teohnidal '.Me.m6 : ranjiufl},v.Nov' 1005 



The . integration gives . 

Mm = Pi* hi J i-Csin-'yCe sin.cp) 



• bh 
+ P 2 -g 



- i j^cos cp J\ — G 



3 sin 2 cp +' | (1 - e 2 ) sinh 1 



e cos cp 



/l .- G 3 



r -\ h / ' — *~ " 

+ p 3 b h 2 j^i cos cp + j_ sin cpj- - p 4 b k -- J 1 - e 3 sin 3 cp (5.3) 



and, with (3.7) added, 



M = 
-S 2 



h i (cos 



cp J\ — 



e 2 sin 2 cp 



1 - € 



£ cos cp 



+ h 2 e cos cp 



+ sinh 

€ 7l - e 3 

+ i hirCsin -1 ) (e sin cp) + eh 2 sin cp 



+• 2 k h~/l - 



e 3 s in 2 cp> 



(5U) 



SJince these moments are small' compared to (3.9) in the ratio 
hj h 2 

— and — , a rough estimate in which only the lowest 

b b , . . .■ • 

powers of e are retained, is sufficient. 

Acc or dingly , 
pab 



M 



m 



+ e h 2 ) cos 



cp + j (hi + e h 3 ) + 2 k h s J- 



or 



pab 



M T (_m) = — — - [(-hj' + 2'g h 2 ) sin cp cos cp] 

M^m) = - [h x cos 2 - cp + (h x + eh 3 -)] 

2 

pab _ 
^s(SL) - — =-* [- ehj sin^ cp cos 3 cp + 2 h 2 ] 
2 



(5.5) 



Since at a > b the component M 2 is greater than. 
Mj , " according to figure 11, the inertia moment in M 2 
direction itself will be enlarged on the ring dimensions; 
hence it is logica.l to assume that h x a.eh 2 (fig. 13). 
With the subsequent simplifications 
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' Mm (m ) = ^ESiS. c "h a s i n cp cos cp 

2 , 

M j (m ) =. . --- ' C h g c o s 8 cp 
3 

M 2 (m) = pati ti s 



(5.6) 



,i.t. is readily .'seen r t : hat" th'ese quantities'' caSi' have no effect 
on the determination X, since M : ' had acted no part 
previously and M- and M s "both are-., small in the ratio 
€ a -h. a ■ 

compared to the previous values. The maximum amounts 

"b 

themselves .change, .in a.ll- three components by quantities 
that are small in the ratio- h 2 /b — hence inclusion of the 
moments due- to- eccentr-ic load application are not worth 
while ', . which leaves for the 'ring dimensions the extreme 
yalues of. figure 12 ,, supplemented • perhaps by a small safety 
margin. ; - • • 



LOADING- -IN 'PURE. TENSION a 1 i 



The calculating process is. the same as before. In the 
axial— force .table (3.6) p x - = a 1 t and p 2 = P3 = ft = 0 > 
for the moment, according to (3.8): 

M ..= .a r t b 2 / j (sin -1 ) (esin cp) [j sin ^ - k cos \|/cos cp] —--^ 

(6.1) 

The .integrations again give 



M„ = a j t b : 



Jl — , -e 2 sin 4 cp Vy fc 



3 ( -r (s in 1 ) (e sin cp) 



.+ X sin I 



3 1 I 



M l - G i * 



jsin s cp - s'in cp <A - e 3 's in 2 'cp ^ i-( s in 1 )"(esin cp^ 



M 2 = - a 



1 t ^\ l-S 

. E3 



2 . 4 

e sin cp 



+ sin 3 cp cos \J/ : ^|-'( c as r ) (e s in cp)^j ' + X < 
or, developed for small 'e : 



cos \|/ 

(6.1 ' ) 
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M m = a 



M, « 



M a = 



jtb 3 j j^sin cp. cos cp ^1 



+ %- sin 3 cp + %~ sin cp + . 
6 2 



— £ X s in cp cos cp r 



a 1 tb"* e * 



a lt D» ^ 



3 
r 



15 



sin 



^1 + ^ sin 4 cp +. . ^ 
cp + . . . J 



3£ 4 



£ L 2 
+ X 



0-T 



sin e! cp cos ' 



cp - . . . j 



|^1 +. ~ s in 4 cp — ~ sin 6 cp + ~— s in 8 cp+ . . . 



(6.1») 



The maximum moments are obtained according to the pre- 
vious considerations. ifhey are 'shown in figure 14. In the 
vicinity, of £ =0, that is, on a ring that does not depart 
too much from the flat ring, the curves (figs. 12 and 14) 
are in complete agreement — the loading is, indeed, approx- 
imately the same except for the 90° rotation of load direc- 
tion. In the region £ ~ 1, on the other hand,. the stress 
of the ring is typically different in the two load cases: 
at point cp = 0 the ring is smooth, hut at cp = tt/2 it, 
has a distinct precurvature which at £ = 1 degenerates 
into a discontinuity; hence the load direction (and, of 
course, the subsidiary effect of the small cylinder in the 
first load case) is essential for the type of stress in 
point . 

Possible approximating formulas, which fail, however, 
in this instance near £ - 1, are: 



M. 



M. 



Mr, 



'max 



max 



max 



4 1 

e »f 2 

- CTitb 2 1 + - £ 

3 1 V 5 



12 1 



(6.2) 



LOADING OF LARGE CYLINDER IN PURE SHEAR 



With T. 



= T 



= T 



- • - indicating the two conjugate 
stresses, the first is applied at an area x = const.; 
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S3-. 



(-hence Wi thr- normal: • and' fails in" ve'eVor.' 'direct i bh' 



1 — sin 2 cp — e £ sin- cp). 

hence stres'ses "ihe ring-, aceop'dlng' to-(2.8), with a load 
per unit length of •. 



cos cp . 

Pi = y=---= ======= (j VI - e^sin 2 cp - c k sin cp) (7.1) 

v 1 - e 2 s.in. . cp . - ^ . - 

The second stress is applied at an area with the normal 
t and falls in direction i ; hence, according to (2 . 8 1 ) , 
Stresses the ring by ~ 

£s = L sin -9-C!0 s ^. • . (7,2) 

The total loading of the ring is p = p : + p 2 , whence, 
according to (3.2) follows the "shear load IT in the form 

I =T tb^ i cos cp - j sin cp - j§ -J\ -.£2 g -in 2 cp y + £ Q (7.3) 

The load distribution being-ant i symmetrical , this 'time .with 
respect to cp = 0 and H- , the axial load N must disap- 

pear at points 0 and 2 for reasons of symmetry; whence 

2 

" • - S 0 " = J X k (7.4) 

b 

Quantity X remains statically indeterminate. From (3.8) 
and (7.3) is obtained ' 



-t X H = Ttb i [*i 

sin" cp 



6 [- JTZprrrr 



. / . , r g r~s — e 2 s in cp c os 2 cp \ 

+ _j sin cp cos V VI in cp + - : : =k cos>|/j 

V v 1 — e 2 s i n 3 cd ' 



e "^s in^cp 



+ k(cos 2 cp — sin 2 cp ) c o s \\i - '- X co's^ (i cos cp + j sin cp) 

b 

and then, according to (3.3)', 



M = T tb 2 j i -- (s in -1 ) ( 6" s in cp) - j c os • cp i "J 1 - e 2 s in 2 cp 
1_~ € 2 • - ' € 

+ k(sin cp cos cp) - X (i sin cp - j cos cp) + M Q "| (7.5) 
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The integration, constant M Q is z.ero -s i'nce. the two bend- 



which must "be antisymmetr i- 



ing moments and M^ o) 

cal with respect to' > = 0 and - disappear. With (2.5) 
and (7.5) -the three, components", of the ..moment .are : 

lTl 



. . s 

cos cp "] 

+ -==r=™~=~ + X cos ^ 
«/l - e 2 s in cp _ 

Mj = - Ttb 2 - sin cp(sin -1 ) ( c sin cp ) 

-sin cpcoscpvl- 



€ 2 sin 2 cp 



M s = | tb 2 — ^ sin cp sin.ty (sin 1 ) ( e sin cp) 

* C € 1 

+ sin 3 cp cos .cp cos\|/j- X sin V 

_r 

or developed again' for small £ : 



M 



(o) 



Ttb 2 i j 1 + sin 4 cp 



4/5 a 
+ G ! — sin cp 



M 



M 



(0) 
1 

(0) 



Ttb 2 € 



2 t 
S + 5 



24 
2 



15 



sin cp 



.in 2 cp 



sin cp cos cp 



Ttb e sin cp cos cp — + — - sm rt cp 

L 3 5 



+ — sin cp + ... 

3 



It is observed that, the statically-, indeterminate 
(which this time is a cross force rather than a foment) 

.nc. M^°) i.s..n,o longer 



does not reappear in 
small, with relation t 



M 



-1 ; whereas 
' (0) 



M. 



so that in the execution 



of the statically indeterminate calculation the. effect 
of M a cancels only for the case of non-oblique bending. 
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When restricted to this particular case (J 13 = 0), the 
calculating' process concerning the determination of "X 
also remains the same as "before, with the sole difference 



that this time the M 



(o) 



portion contains 'the lower e 



•'-powers,' so that the extreme case 



G-Jm 

a = — ± -ss. 0 
3J 2 . 



also adjoins the case a ^ 0 without "discontinuity. 

The formula for predicting X can again he written 



with 



X = 



k x + B 2 a 



0 j_ + D a a 



TT/c' 



TT /i 



A. i = 



M 



.°) 



d cp, B 



tan + d cp 



> (7.6) 



o 
tt/e 



tt/2 TT 

= / cos\l l dcp,n i =/ 



sin 3 \[/ 
cos \j/ 



d cp 



J 



for e > l/2 the integral must he again numerically eval- 
uated; for small e ■ 

2 



X 



l + LI + £ 



1 + IQco c 4 



16 



64 



384 



1 _ £ 0.0175 e 4 +--- + 0.4063 ---- 

16 8 8 



1 - a 



.3 



e 

100 



+ i_-_i±' e 2 + (4.09 -i 4.81 a + 1.56 a 3 ) +..} (7 

8 ' 



Y 

7 



The result is shown in figure 15. This time 

max 

falls in the symmetry points cp = 0 , tt/2, M 2 near the 

mo, x 

is seen to he little 



M. 



point cp = tt/4. Then, 
greater than in the extreme case e 



-max 



at a = 0 ; 



whereas 



M n 



increases considerably near £ = 1. 



The 



'max 
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torsion moment has th«" opposite sign at 0 and tt/2 . At • 

cp = tt/2 the dependence of the maximum torsion moment on 

a is -as expected; it rises with increasing torsional 

stiffness. This aspect of M T (a) is due to the fact 

that the statically indeterminate port ion „. which at cp =. 

tt/2 is smaller- in amount than the statically determinate, 

decreases with increasing a. — At cp = 0 the conditions 

are reversed: Since the statically indeterminate portion 

governs the sign in this instance, (Mm) decreases 

x extra 

with increasing a. The curve of the other "bending mo- 
ment M x has "been omitted in figure 15, since it. is al- 
most straight and would intrude, moreover,, in the range 
covered "by the M T curves. 

The following simple approximate formulas remain: 



1 3 

M 2 = - t th 3 
s max 2 



M, = - Tth s (l + 0.15 c s ) 
max 6 



Mrp 



= _ Tth 2 (l + 0.3 e 2 ) 
'max 6 



(7.8) 



the last one fails near G = 1. 



•Translation "by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Jigs. 1,2,3,4.5,6,7,8,9 



Figure 1.- Section of two cylinders. 
4K 





ligure 2.- The first load case 
, within the limit 

1 «s € >0. 

a 



£t*d£{ 



Figure 3.- Identification of axes 
t, n, It on the flat 

ring. 





- — * 



A/ 



Figure 4.- Equilibrium of ring 
element d s. 



tttt tftttttttttt 



Figure 5.- Equilibrium of flat 
ring element. 





Figure 6.- Resolution of skin stresB 

applied at ring 
element d s.. 



Figure 7.- The three projections 
of the space curve. 





Figure 8.- Defininition of angle ty. 



K 

Figure 9.- Solution of linear load 

£ ( the plot lies in 
the tangential plane of the 
cylinder. 
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Tigs. 10,11,12,13 



til 





pb'a-QS 



9 



a) 4* 

Figure- 10.- Relative position of 

the different normals n. 



pb'a-M 





Figure 12.- Maximum moments 

against h 
parameter GJif € = a » 

Load: internal pressure p. 



IP W «' »* K' T Sf 

Figure 11.- Moment distribution 

along op. The statically 
indeterminate portion X plotted 
against u m Jj£r . (axes ratio e m h m 1 ) 
i f SJ2 a 2 

( internal pressure p.) 



Figure 13.- Bing section, 
moment lever 
arms h^ f h2> 




figure 14.- Maximum moment b against 
D parameter GJqi 
c - 5> « = TSJ Z 

Load: tension Ojt.On comparison 
with figure 12. -note that in the 
extreme case e— >0 the loads v a 
and o\t mutually agree - 2 




figure 15.- Maximum moments against 
d parameter QJrf 
e = a' * " K2 

Load: shear t t in large cylinder* 
M 2 (*r /4)*M 2 ^ increases hy 

decreasing <*. 



